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Abstract 

We find a general expression for the free energy of G(<j>) = 4> 2k generalized 2D Yang- 
Mills theories in the strong {A > A c ) region at large N . We also show that in this 
region, the density function of Young tableau of these models is a three-cut problem. 
In the specific </> 6 model, we show that the theory has a third order phase transition, 
like <ft 2 (YM2) and 4 models. We have some comments for k > 4 cases. At the end, 
we study the phase structure of (ft 2 + gcj)^ model for g < ^ region. 
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1 Introduction 



The pure 2D Yang-Mills theory (YM2) is defined by the Lagrangian tr(F 2 ) on a compact 
Riemann surface. In an equivalent formulation of this theory, one can use itr(BF) + tr(5 2 ) 
as the Lagrangian, where B is an auxiliary pseudo-scalar field in the adjoint representation 
of the gauge group. Path integration over field B leaves an effective Lagrangian of the form 
tr(F 2 ). 

Pure YM 2 theory on a compact Riemann surface is characterized by its invariance under 
area preserving diffeomorphism and by the fact that there are no propagating degrees of 
freedom. Now it is easy to realize that these properties are not unique to itr(BF) + tr(B 2 ) 
theory, but rather are shared by a wide class of theories, called the generalized Yang-Mills 
theories (gYM 2 ). These theories are defined by replacing the tr(B 2 ) term by an arbitrary 
class function A(B) [10]. Besides the reasons which was discussed in [1] to justify the study 
of gYM 2 , we can add two other reasons. The first one is that the Wilson loop vacuum 
expectation value of gYM 2 obeys the famous area law behaviour, [11], and this behaviour 
is a signal of confinement which is one of the most important unsolved problem of QCD. 
Second, the existance of the third order phase transition in some of the gYM 2 theories (one 
case was studied in [5] and other examples will be studied in this paper), is another indication 
for equivalence of YM 2 and gYM 2 as a 2D counterpart of the theory of strong interaction. 

In ref.[10], the partition function of gYM 2 has been calculated by regarding the gener- 
alized Yang-Mills action as perturbation of topological theory at zero area. In ref.[l], the 
partition function has been obtained by generalization the Migdal's suggestion about the 
local factor of plaquettes, and has been shown that this generalizatin satisfies the necessary 
requirements. And in ref. [11], this function has been calculated in continuum approach and 
by standard path integral method. The gYM 2 theories can be further coupled to fermions, 
thus obtaining QCD 2 and generalized QCD 2 theories [1]. Recently there has been a lot of 
interest in these theories. Phase structure, string interpretation and algebraic aspects of 
these theories are reviewed in [2]. 

In this article we explore the issue of phase structure of the gYM 2 theories. An early 
study of phase transition of YM 2 in the large- N limit on the lattice [3] reveals a third order 
phase transition. The study of pure continuum YM 2 in the large- N limit on the sphere also 
yields a similar result [4] . This result is obtained by calculating the free energy as a function 
of the area of the sphere (A) and distinguishing between the small and large area behaviour 
of this function. In [5], the authors considered gYM 2 in the large- N limit on the sphere and 
they found an exact expression for an arbitrary gYM 2 theory in the weak (A < A c ) region, 
where A c is the critical area. They also have found a third order phase transition for the 
specific 4 model. 

In this paper we discuss the issue of phase transition for a wider class of theories. In 
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section 2 we review the derivation of the free energy and the density function in the weak 
(A < A c ) region. In section 3 we study the 4> 2k theories. First, we show that the density 
function of these models has two maxima in the (A < A c ) region, similar to 4 model. This, 
then enable us to use the method of ref. [5] to obtain a general expression for the free energy 
for these theories in the strong (A > A c ) region. In section 4 we compute the free energy 
near the transition point for the specific 6 model and we show this model also has a third 
order pahse transition. We also make some remarks for k > 4 models. And finally in section 
5 we study another class of models, namely the 2 + g</> 4 models. If g < j then the density 
function will have only one maximum at the origin. We show that these models also undergo 
a third order phase transition in this domain. 



2 Large- N behaviour of gYM 2 at A < A c 

The partition function of the gYM 2 on a sphere is [5, and references there- in] 

Z = Y.dle- AK ^\ (1) 

r 

where r's label the irreducible representations of the gauge group, d r is the dimension of the 
r-th representation, A is the area of the sphere and A(r) is : 

^■) = t^hC k {rl (2) 

k=l 

in which Ck is the A;-th Casimir of group, and a^'s are arbitrary constants. Now consider the 
gauge group U(N) and parametrize its representation by n\ > n 2 > ■ • • > n N , where rij is 
the length of the i-th row of the Young tableau. It is found that 

d r = n 

l<i<j<N J 1 

N 

C k = Y^[(ni + N-i) k -(N-i) k ]. (3) 

i=l 

To make the partition function (1) convergent, it is necessary that p in eq.(2) be even and 
a p > 0. 

Now, follwing [4], one can write the partition function (1) at large N, as a path integral 
over continuous parameters. We introduce the continuous function: 

(f)(x) = -n(x) - 1 + x, (4) 

where 

< x := i/N < 1 and n(x) := m/N. (5) 



2 



The partition function (1) then becomes 

Z= [ J] d(f)(x)e SWx)] , (6) 

where 

S(0)=N 2 {-A f 1 dxG[</>(x)] + C dx f'dy log\<f>(x) - 4>{y)\ }, (7) 
Jo jo Jo 

apart from an unimportant constant, and 

G[<f>] = E(-l) fc a fc fc . (8) 
k=i 

Now we introduce the density 

where in the cases which G is an even function, the normalization condition for p is 

P p(X)d\ = I. (10) 

J— a 

The function p(z) in this case is [5] 



*> = ^| o 2 4r + 1 , ) ; 1)! ^^-'(°). (id 

where 

m = ^g' (</>), (12) 

and 

g(n) is 

the n-th derivative of g. Similarly one can express the normalization condition, 

eq.(10), as 

t^^^m^. d3) 

Defining the free energy as 

F := ~\nZ, (14) 
then the derivative of this free energy with respect to the area of the sphere is 

F'{A)= f 1 dx G[<f>(x)]= r dX G(X)p(X). (15) 

JO J -a 

The condition n\ > n 2 > • • • > imposes the following condition on the density p(A) 

p{X) < 1. (16) 

Thus, first we determine the parameter a from eq. (13), then we calculate F'(A) from eq. 
(15). Note that the above solution is valid in the weak (A < A c ) region, where A c is the 
critical area. If A > A c , then the condition p < 1 is violated. 
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3 The G{<t>) = (f) 2k model 

In order to study the behaviour of any model in the strong (A > A c ) region, we need to 
know the explicit form of the density p in the weak (A < A c ) region. From eq. (11) we can 
obtain p for any even function G{4>). However, in this section we consider a simple case, 
namely the G((f>) = <p 2k model with arbitrary positive integer k . The density function p in 
the weak region is 

The interesting point is that the above density function has only one minimum at z — 0, 
and two maxima which are symmetric with respect to the origin. To see this, notice that 
equating the derivative of p equal to zero, will yield z — 0, and 

f(y) = -i + E ( !?/~ 1)! ;, y~ (n+1) = o, (is) 

where y = z 2 /a 2 . Now as all of the coefficients of y, in eq. (18), are positive, the function 
f(y) is a monotonically decreasing function and has only one root,ie. y . Next, expanding 
p near the origin, we obtain 

thus, p"(0) > 0. The origin is then a minimum. But, near the points z = ±a the curve 
p(z) is concave downward, consequently the points z^ = ±a^/y^ will correspond to two 
symmetric maxima of the density. Therefore in the strong region, all 2fc 's models are three- 
cut problems. The function F'(A) for G(4>) = <f) 2h in the weak regoin is [5] 

FUA) = jjL. (20) 

Next, we study the strong (A > A c ) region. Following [5], we use the following ansatz 
for p 

p.(z) = l 1, ^[-b,-c]\J[c,b]=:L' 

\ p s (z), z G [-a, -b] U[-c, c] \J[b, a] =: L. 

Then, if we define the function H s (z) [6] 

H s (z) := P f dw (22) 

J —a Z — W 

where, P indicates the principal value of the integral, it has the follwing expansion at large 
z 

H s {z) = - + - 3 P Ps (\)\ 2 d\ + ■■■ + -^Fi(A) + ■■■■ (23) 

Z Z J —a Z 
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Hence, one can obtain F' S (A) via expansion of H a (z). 

To calculate the function H s (z), we repeat the same steps which was followed in [5], and 
using some complex analysis techniques ([7]), which yield the following result for <p 2k model: 

00 r b \rl\ 

H s {z) = kAz 2k ~ l — kAR(z) £ a(n,p,q)z 2 ^ n - p -^-2R(z) - (24) 

where 
and 



= ^/(z 2 -a 2 )(z 2 -b 2 )(z 2 -c 2 ). (26) 

Now, we expand H s (z)/R(z) and demand that it behaves like l/z 4 at large z. It can be 
shown that the coefficients of positive powers of z in the above expansion are identically 
zero. Next, we calculate the coefficients of l/z 2 this gives 

00 r b Ar/A 

kA £ a(n,p,q)-2 — = 0, (27) 

n,p,g=0 Jc ^H A J 

in which n + p + g = /c — 1. By setting the coefficient of l/z 4 to unity we obtain 

M £ a(n,p,q)-2 ^ = 1, (28) 

where n + p + q = k. In k = 2 case, 4 theory, eqs.(27) and (28) reduces to 

- 4 < a2+62+ ^ 2 /H' (29) 

and 



-(a 4 + 6 4 + c 4 ) + -(a 2 6 2 + 6 2 c 2 + c 2 a 2 ) 
4 2 



which are the same equations that was obtained in [5]. 

We can express the action in terms of p s (z). If we maximize this action along with the 
normalization condition eq.(10), as a constraint, we obtain another equation. The procedure 
is fully explained in [5,8]. The result is 

kA E f "(n,P, q)z 2{t -'-'-- 2 '>R(z)dz + 2 fdz pf J^L ^ W 



n,p,9=0 



Note that we can determine three unknown parameters a, b, and c from equations (27), (28), 
and (31). 
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To compute the function F' S (A), we start from the function H s (z) directly. First we 
expand the R(z) 



R(z)=z 3 d(l-- 2 )(l- U -)(l-- 2 ) = -z" Y, P{n',p',q')z~ 2{n ' +P ' +q '\ (32) 

V Z Z Z I I I n 

' n,p,q'=0 

where 

= ( 2n-Zn2,-Z)m-Z)« 
2 n +P + in\p\q\ 

Furthermore, we define (—3)!! = — 1. If we substitute the above expansion into the eq.(24), 
we obtain 

H (z) - kAz 2 ^ 1 + kA V a(n,P,gM<p / ,g / ) 

JTL S \Z) — hJ±Z f hJ± ^ z 2( n+p+q -k)+2(n'+p'+q')+l 

n,p,q,n',p',q'=0 

From eq.(23) we see that the coefficient of l/z 2k+1 in the expansion of H s (z) is F' S (A). 
Therefore from the above equation we get 

F' s (A) = kA £ a(n,p,q)(3(n',p',q')+2 £ (3(n',p',q') J ^— ^. (35) 

n,p,q,n' ,p' ,q'=0 n,n',p',q'=0 c \ / 

In addition, in the first summation of eq. (35) we have the follwing conditions on the indices: 

(n + p + <?) + (ra' + p' + q) = 2k, (36 - a) 

and 

2A; - 2n - 2p - 2q - 4 > 0. (36 - b) 

The condition (36-a) appears due to selection of a specific power of z in the expansion, and 
the condition (36-b) is due to complex integration. Furthermore in the second summation 
we have the follwing condition on the indices 

n + (n'+p' + q') = k + 1. (36 - c) 

In this way, we find the explicit relation of the free energy of 2fc 's models. For the k = 2 
case, our results agree with those of ref. [5]. 
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4 The G(</>) = (f model 

As an application of the previous results, we will investigate carefully the phase structure 
of the G{4>) = 6 model. From eq.(17) we obtain the density p for this model in the weak 
region, the result is 

P(z) = — (-r- + ^T + z*)^^- (37) 

71 O I 
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From the normalization condition ( eq.(13)) we obtain a = (16/ '(15 A)) 1 / 6 . In addition from 
eq.(20) F^(A) = 1/(6 A). This density function has a minimum at z=0 and two maxima 
at = ±(\J + l)a/2. At A = A c , the density function (37) is equal to unity at . 
Using this fact, we find the critical area A c 

3125 156W3 
Ac = 7V W-^3312")- (38) 

For A > A c the eq.(37) is not valid. 

Next, we analyse this model in the strong (A > A c ) region. The equations (27), (28), 
and (31) in this case become 



3A 



h a * + b 4 + c*) + ](a 2 b 2 + b 2 c 2 + c 2 a 2 ) 
8 4 



: b \d\ , . 

- 2 /fl(A)= ' < 39 > 



3A 



and 



A( a 6 + 6 6 + c 6) + 3 ^4 + + + ^ + ^4 + ^ + 1 w 

16 16 8 

To study the structure of the phase transition, we must consider the theory on the sphere 
with A = A c + e area, where e is an infinitesimal positive number. In this region, following 
[5], we use the follwing change of variables 

c= s(l-y), 

b = s(l + y), (42) 



a = 5^2^-2 + 



Note that these parameters are introduced so that at critical point, e and y are equal to zero 
and s is equal to Zq. Now expanding the equations (39), (40) and (41), we find 

Ml8 _aV3)-iL + (M^-3) + J(i + |))e 

+ (^-^ + f))«' + ((9 + ^-^(| + ^)y 

3 As 4 u 11^3 89 xx 9 /n „ 4 ?r ,5363 129^,, , x 
+ '— + + + (3 - 4s " ^liH + — W = °" (43) 

and 
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(39v/3 - 57)As e - 1 - - + ((42 - 18V3)As e + -(]- + %)e + ((^5 _ 

T] 2 6 8 2 

+ ^ + ((33 + 3V3)^ 8 - ^(2v/3 + ^)), 2 + ((? + ^)As* 
7rs.35V3 121,. 2 ni . . , 7ts. 199^3 8267,, , n 

+ - (-T" + x ))ey + ((3v ^ + 24)As ~ + T92" ))2/ = °' (44) 



and 



3(l + V3)^-^ + f) 



(( 



5v/3 



j i 2y3 ,,7v3 13, . = 



7T.5 , 2 , ,n , a /o\ ^ r vr^lv^ , 211,, 2 , „25 , 15^3 



+-(- + ^))e 2 -((2 + 4 v / 3>L S 5 + . 
7] 2 9 r] 6b 



+ 48 )y + (( T + 2 



71,635 275\/3 

+ -( 1 — 

if 24 18 



))e|/ 2 - ((- + 



56 32^, tt. 60673^ 23353 



-) + -( 



+ 



■))y A = o. (45) 



3 3 ' r] y 1728 384 

The parameter 77 = \]2y/% — 3 is used for the sake of brevity in the above formulas. We also 
have kept terms up to order y 4 or e 2 (we will show that e is of order y 2 ). Next we obtain s 
from eq.(43), the result is 



.7T.1 .12 + 7^,^ 
S = ( — 5 ( ) 10 

K A J 1 648 ; 
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'32 ' 1~2 



A V3 X 2 ,167 95^, 2 ,911 llVS, 4 

+(7 + - (ttt + ^7-)^ + (7777 + — r-)f 



^4 ' 6 '96 
Substituting this s in eq.(45) will give 



96 



192 



(46) 



( 



5v/3 1, 2 ,15 37v/3, 4 
- ' )y . 



16 + 48 



(47) 



So e is of order y 2 . Using eq.(44), we obtain 



AV3 2 317 77yg 2 
v 15 5 } ^360 150 ' ' 



(48) 



and 



,11 17^ 

e = ( y-^5- )5 



^533^ 14929^ 2 



The parameter 5 is the reduced 

-A, 

From eq. (35), we find F' S (A), the result is 



3600 

A-A c 



3600 



(49) 
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F' S (A) 



3A 



35(a 12 + b 12 + c 12 ) - 12(a 6 6 6 + a 6 c 6 + 6 6 c 6 ) 



1024 

+12aW(aV + a 4 c 2 + b 4 a 2 + b 4 c 2 + c 4 a 2 + c 4 b 2 ) + 14aVc 4 
-2a 2 b 2 c 2 (a 6 + b 6 + c 6 ) - 19(a 8 6 4 + a 8 c 4 + b 8 a 4 + b 8 c 4 + c 8 a 4 + c 8 b 4 ) 
-10(a 1( V + a 10 c 2 + b w a 2 + b 10 c 2 + c w a 2 + c w b 2 ) 



+ 



A( a 8 + & 8 + c 8) _ l( a 6 6 2 + fl 6 c 2 + ^2 + ^2 + ^2 + ^2) 

64 16 

+ l a 2 6 2 c 2 (a 2 + ^2 + c 2 ) _ 1 (fl 4 6 4 + ^4 + ^4) 

1 o o.z 



1(A) 



+ 



8 



a 6 + 6 6 + c B - (a 2 6 4 + a 2 c 4 + 6 2 a 4 + 6 2 c 4 + cV + c 2 6 4 ) + 2a 2 6 2 c 



2i4 i „2„4 i r,2„4 , 7,2 „4 , „2 „4 , „2l4\ 



,2i,2„2 



b X 3 dX 



R(X) 



+ 



a* + tf + c l -2(a 2 b 2 + a 2 c2 + b 2 c^ 



b X 5 dX 



+ (a 2 + b 2 + c 2 



[ b ^dX_ [< 
Jc R(X) Jc 



X 9 dX 



R(X) v ' Jc R(X) Jc R(X)' 

(50) 

To compute F^(A), we express the parameters a, b and c in terms of 5, and after a lengthy 
calculations eq.(50) reduces to 



. 271 1289^3, r2 



10800 504000 

If we compare this with F' W (A) given in the begining of this section we find 



(51) 



F' S {A) - F' W {A) 



271 1289^, 1 ,A-A C1 



'64800 3024000 M r A r 



Y + 



(52) 



Therefore, we have a third order phase transition, which is the same as the ordinary YM2 [4] 
and G(<f)) = <fi 4 model [5]. In principle one could study other <p 2k models in the same manner. 
The k = 4 case is solvable analytically, but the expressions become too complicated. For 
k > 5 one can not solve the equations analytically. 



5 The G(4>) = (j) 2 + g<j> A model 

So far in our study of the phase transition for gYM 2 theories we considered only those 
G(4>) which contained a single term. In ref. [9] the authors study the phase transition of 
gYM.2 on a closed surface of arbitrary genus and area A. In particular they investigate the 
G(4>) = (p 2 + g<p 3 model. However, their treatment is mostly qualitative. In this section we 
consider a simple combination of <fi 2 and 4 , namely we study the G{<p) = <fi 2 + g<fi 4 model. 
In the weak region we can obtain the density p from eq.(10), the result is 

p(z) = -Va 2 - z 2 {l + ga 2 + 2gz 2 ). (53) 
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The above density will have only one maxima at z — 0, when 

3ga 2 < 1. (54) 

The normalization condition (13) yields 

X -Aa 2 + ^gAa A = 1. (55) 
Using (55) the condition (54) reduces to 

g < A/A. (56) 

Therefore, if this condition is satisfied, we will have a two-cut problem in A > A c areas and 
we limit ourselves to this region ( condition (56)). 
The free energy of this model is ( using (15)) 

K(A) = \a'A + ^ga«A + ^VA (57) 
To study this model in the strong (A > A c ) region we use the following ansatz for p [4]. 

[Ps{z), z E [-a,-b\\J[b,a\ 

Using complex analysis [4,5,6] we obtain the function H s (z) defined through eq.(22), the 
result is 



H s (z) = Az + 2gAz 3 - ^/(z 2 - a 2 )(z 2 - b 2 ) 2gAz + J 



b (z-X)U(X) 



(59) 



where 



U(X) — sj{a? - A 2 ) (ft 2 - A 2 ). (60) 

Using eq.(22) it is seen that H s {z) should behave as \ for large z. Therefore the coefficient 
of z of eq.(59) must be equal to zero 

A+ ° AM -lm =0 - (61) 

and the coefficient of 1/z must be one 

-MA + gA(-M 2 -N)- t ^ = I. (62) 
2 y V 4 ; J-bU(X) V ; 

In these relations M = a 2 + b 2 and iV = a 2 b 2 . The two unknown parameters a and b can be 
determined from these two equations. 
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Using equations (15), (22) and (59) we obtain the following expression for the free energy: 
= 8 ~ 2^ L m + 2 M Lu(X)- L 



+9 



b U(X) 

'--MN + —M 3 )I T77TT + {-M 2 --N) T77TT + -M 

8 2 J — 



^ A 4 dA 




-6 C/(A) 




/• b A 4 dA 


A 6 dA" 


7-6 Z7(A) 


7-6 Z7(A)_ 



4 16 J J-bU{X) v 8 2 y 7- b f/(A) 2 

- 2gA(±MN - ^M 3 ) - 2g 2 A(~N 2 + ^M 2 A - ^M 4 ). (63) 
In order to investigate the phase transition, we use the following change of variables 

a = a c (l - h), (64) 

(-) 2 = k. (65) 
a 

The parameter a c is the value of a at the critical point, and h and k are equal to zero at this 
point. At the critical point, the normalization condition (55) becomes 

X -A c a 2 + -gA c a c 4 = 1. (66) 

In addition, as at A = A c , p(0) is equal to one, we obtain 

a c A c + gA c a c 3 = ir. (67) 

In order to study the phase transition of this model, we limit ourselves to the small values 
of g (g « 1). In this limit A c and a c become (by eqs.(66) and (67)) 



7T 2 Aq 2 

A c = — -g+^ + ---, (68) 



and 



a c = + — — + • • • . (69) 

Now, we follow the same steps of previous section. First, we expand the equations (61) and 
(62) in terms of h and k and then solve these equations for h and k. The results, up to order 
5 2 , are 

„ (2 + 14 A + 42 J, + ( _Z_|A + -^ + .... (71) 

Second, we find F^(A) from eq.(63). The result is 

^) = 2 {- 2 - -A + 1% + d + 3 f + 85 ^ 2 + • ■ •>■ < 72 > 
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It is understood that this relation is true only up to order g 2 and S 2 . To compare the equation 
(72) with F'(A) in the weak region, we first note that a 2 , up to order g 2 , is (from (55)) 

a2 = l( 1 - 3 f + 18 f^)' ( 73 ) 
and from this, F^{A) becomes (from (57)) 

K(A) = \(\-\^ + %). (74) 

Therefore 

F' S {A) - F' W {A) = + sf- + 85^)(^-^) 2 + • • • . (75) 

This shows that the G((f>) = <p 2 + g0 4 model in g < A/A region has also a third order phase 
transition. Also note that at g — 0, the eq.(75) reduces to 

F1(A) - FUA) = - 2 (^-^r + ■ ■ ■ , (76) 

which is the same relation that was obtained in ref. [4] for ordinary 2-dimensional Yang-Mills 
theory. 

At the end we should mention that for g > Aj '4, the density function p (eq.(53)) has two 
symmetric maxima in the weak (A < A c ) region. So we have a three-cut problem and the 
above perturbative calculation is not possible. 
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